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i Abstract 

Using the coherent state techniques developed for the analysis of the EPRL model we give 
the asymptotic formula for the Ponzano-Regge model amplitude for non-tardis triangulations 
of handlebodies in the limit of large boundary spins. The formula produces a sum over all 
O , possible immersions of the boundary triangulation and its value is given by the cosine of the 

Regge action evaluated on these. Furthermore the asymptotic scaling registers the existence 



5-H , of flexible immersions. We verify numerically that this formula approximates the 6j-symbol 



bJO, 



o 



for large spins. 



> ■ 1 Introduction 

In [1], Ponzano and Regge gave a formula for the large spin limit of the 6j symbol. It was 
^SJ ' found to be related to the Regge action for discrete general relativity and with this motivation 

. they constructed the hrst spin foam model of 3d gravity. Their asymptotic formula was first 

I proved in [2] then more recently using different methods in [3l Hj and the square of the 6j symbol 

was also studied in the context of relativistic spin networks [3 [6] . The next to leading order 
approximation was recently considered in . The precise formulation of the full state sum was 
studied in |8j. 

^ ' In [HI dn] , the semiclassical limit of some recent spin foam models WI[ [T^ was analysed 

' using the coherent state techniques introduced in [13]. In particular the boundary there was 

formulated in term of coherent tetrahedra. Here we apply the same techniques in the 3d case 
using coherent triangles and, instead of a single vertex amplitude, we analyze triangulations of 
arbitrary genus handlebodies. This finally opens up the possibility of a continuum limit and 
renormalization analysis of the model for this restricted class of 3-manifolds. In particular the 
resulting formula is well suited for studying the graviton propagator as introduced for Ponzano- 
Regge in [Ti] . 

We begin the paper by describing the formulation of the Ponzano-Regge model in terms of a 
single spin network diagram dual to the boundary triangulation. We then describe the boundary 
state choice in detail and proceed to give the asymptotic formula in terms of immersions of the 
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boundary triangulation. The asymptotic scaling of the ampHtude has the interesting feature 
that it registers whether or not there are flexible immersions of the boundary. An explicit 
example is provided in Steffen's polyhedron. This analysis sheds light on the general way in 
which asymptotic behaviour for larger triangulations can emerge from a spin foam model. In 
particular it does not need to proceed by taking the asymptotics of the individual simplex 
amplitudes first. 

2 The Ponzano-Regge model in terms of coherent states on the 
boundary 

The Ponzano Regge amplitude was originally defined in terms of 6j symbols with a cutoff regular- 
ization on the interior vertices. More recently, it was shown in that the cutoff regularization 
for sums over representations in some cases disallows a 2-3 Pachner move (the Biedenharn-EUiot 
identity does not hold for a restricted sum over representations.) This meant topological in- 
variance of the partition function can not be proved with Pachner moves in this form. The 
alternative formulation in terms of delta functions and integrals over SU(2) regularized with a 
gauge fixing tree is both finite and invariant under Pachner moves. Another regularization using 
representations of a quantum group is given by the Turaev-Viro model, however it was also shown 
in [8] that the limiting procedure that gives the Ponzano-Regge model is only known to exist 
for so-called non-tardis triangulations - i.e. a triangulation whose edge lengths are restricted to 
a finite range by the boundary edge lengths. In order to avoid discussing regularization, in this 
paper we will restrict to only considering these non-tardis triangulations which are by definition 
finite. Slightly extending the terminology of [S] , we will call a manifold a non-tardis manifold 
if there exists a 'non-tardis' triangulation of S^. 

For a 3-manifold with orientable 2-boundary its boundary state space is then given 
by the possible geometric triangulations of the 2-boundary with half integer edge lengths. The 
amplitude for such a non-tardis manifold is given in terms of a non-tardis triangulation T of 
that extends the boundary triangulation, and some boundary state ^: 

zM^, ^) = E n d-(^-e) n n^Tet). (d 

je e A ^ ' a- 

Here e is an edge, A a triangle and a a tetrahedron of the triangulation of the interior, j are half 
integers labelling the irreps of SU(2). The amplitudes (Theta) and (Tet) are the spin network 
evaluation of the theta graph and the planar tetrahedral spin network respectively. These spin 
networks are the two dimensional duals to the interior triangles A and, respectively, to the 
surface of the tetrahedra a in the triangulation. The labelling of the spin network surface duals 
of the A and a is given by assigning the j associated to each edge to each dual edge that crosses 
it. Finally dim(j) = (—1)'^^ {2j + 1) is the (super)-dimension of the jth SU(2) irrep in graphical 
calculus. 

In the interior the normalisation and phase of the intertwiners cancels. However, at the 
boundary these are arbitrary normalisation for each face. This information is in the boundary 
state ^ which consists of the boundary edge length data and the particular intertwiner chosen 
at each face. 

2.1 Ponzano Regge on the boundary 

In some cases it is possible to reformulate the Ponzano-Regge model defined above as a spin 
network evaluation on the 2-boundary of the manifold. In fact, Ponzano and Regge originally 
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constructed the state sum model such that it agreed with the evaluation of a planar spin net- 
work associated to the boundary of a 3-ball. An algorithm to construct a non-tardis interior 
triangulation given an arbitrary triangulation of the boundary of B^, was given using recoupling 
theory by Moussouris in |15j . This algorithm consists of reducing the boundary spin network 
to a product of 6j symbols (which is always possible for a planar diagram) using the recoupling 
identity and Schur's Lemma and then reconstructing the interior triangulation from these 6j 
symbols. Since the boundary spin network is finite, this procedure gives a manifestly finite 
definition of the partition function. 

In this paper we will extend this result to spin networks on the boundary of handlebodies of 
arbitrary genus. A non-tardis triangulation of a handlebody of genus g can be constructed as 
follows. Start with a triangulation of the boundary of with g distinct pairs of triangles that 
do not share a common vertex. The boundary of the handlebody can be formed by identifying 
these triangles and a non-tardis triangulation of the interior is given by applying the Moussouris 
algorithm. This procedure may result in a degenerate triangulation of the handlebody even if 
the triangulation of the ball is non-degenerate. 

We will begin our analysis by reformulating the amplitude for the 3-ball B^ on a non-tardis 
triangulation as a spin network on the boundary by a "reverse Moussouris algorithm." We then 
describe how this procedure is altered for handlebodies of arbitrary genus. From now on, 
denotes a handlebody. 

Lemma 1. The Ponzano-Regge amplitude for a non-tardis triangulation of B^ can he expressed 
in terms of a single spin network evaluation 



where * is the two dimensional dual of the surface triangulation with each dual edge labelled with 
the SU(2) irrep corresponding to the length of the edge it is dual to, and the spin network is 
evaluated as the planar projection without crossings, with the intertwiner normalisation given by 

Proof. In order to reexpress the 3-ball with a given triangulation T and the amplitude 
Z{^,B^) as the spin network evaluation of its boundary we proceed inductively. Note first that 
a triangulation of B^ given by a single tetrahedron is already of the form we want to put it in: by 
([1]) its amplitude is given exactly by the evaluation of the spin network dual to its boundary 2- 
geometry with an intertwiner normalisation chosen at each surface triangle. This establishes the 
base case. We now need to show that the statement remains true when one glues tetrahedra on 
to a ball amplitude already expressed in this manner, and thus reconstruct arbitrary non-tardis 
triangulations of the 3-ball. To glue we add the necessary face and edge amplitudes for the new 
interior faces and edges with the same the normalisation and phase choice of the intertwiners 
chosen in the boundary state before. These boundary choices will therefore cancel. This is in 
accordance with the observation above that the phase choice and normalisation on the interior 
are left arbitrary. A tetrahedron can be glued onto a ball non-degenerately with one or two 
faces: 

1. If we glue one face of the tetrahedron with one face of the ball we create an inner triangle. 
The PR amplitude of the new ball differs from the old one by a ^-j-hcta) ^ tetrahedral 
net. In the spin network evaluation the vertices of the 3-ball and the tetrahedral amplitude 
corresponding to the glued face, together with the face amplitude, are the normalized 
projector on the invariant subspace of the irreps on the edges. As both amplitudes being 
glued already are invariant we can simply replace them with parallel strands, see Figured! 
This changes the spin network graph being evaluated by changing a vertex to a triangle. 
This is the dual to the change of the surface triangulation, and the resulting amplitude 
still satisfies the lemma. 
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Figure 1: Case 1: Reduction of the PR amplitude for two tetrahedra to the spin network on the 
boundary. 




Figure 2: By replacing the gluing along the 4 faces with a group averaging on edges crossing the 
dashed line (left) we reexpress the T amplitude on the boundary (right). 



2. If we glue two faces of the tetrahedron onto the ball, we create an inner edge and two 
inner faces, the PR amplitude changes by adding a tetrahedral net, two thetas and one 
dimension factor. However, these nets correspond exactly to the 6j symbol for changing 
the ball amplitude from being connected along the dual of the old boundary to the dual 
of the new boundary. 

Note that for any non-tardis triangulation of B"^ we can always build it up from a single tetra- 
hedron by gluing on one or two faces. Furthermore the two operations described above do not 
introduce crossings and respect the planar projection chosen. 

This establishes that one can express the Ponzano Regge amplitude of an arbitrary triangu- 
lation of as a spin network evaluation on its boundary S^. This proves the lemma □. 

Consider next the case of a solid torus x S^, which we call T. Take a disc D C T such 
that B = X {p} £ D"^ X S^. For future purposes, note that it intersects a non-contractible 
loop in T. We can now always move this disc by a homotopy that keeps dO on dT such that its 
boundary is the union of at least three edges of the boundary triangulation. Due to triangulation 
invariance we can then choose a triangulation such that B has no internal vertex. We can then 
cut the Ponzano Regge amplitude along this surface, the resulting space is topologically B"^ and 
we can apply the previous lemma. This yields a ball where two discs on the boundary are glued 
by identifying edges and using the PR face and edge weights. Call n the number of edges that 
make up dlD). As we chose a disc with no internal vertex, the spin network dual to it has to be 
an n — 2 vertex string with one outgoing edge per inner vertex, and two at the end. Together 
with the face amplitudes this is simply the projector onto the invariant subspace of the irreps 
on the circle dO. This projection can then be replaced by a group averaging on the strands dual 
to dB: 

Zpn{^,T)=[ dh{{dT)l) (3) 

JSU(2) 
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where (9T)j^ is the spin network dual to the surface of the torus with h inserted along dual edges 
crossing dD. The diagram is defined by first cutting along this circle, choosing the planar no 
crossing diagram of the graph and then connecting up along the identified edges. See figure [21 
and Appendix |A] for an explicit example. 

We can easily generalize this example to arbitrary genus handlebodies. By definition, a 
handlebody of genus g comes equipped with a set of g standard cuts that reduce the handlebody 
to the 3-ball. We call these cuts Bj with an index i G C where C is a set of labels for the standard 
cuts. For later use, note that it is always possible to define a complete set of generators Cj of 
the homology group such that each q is transversal to the cut Dj and does not intersect 

the other cuts. We can choose an equivalent set of cuts that are related to the standard cuts by 
boundary preserving homotopy as long as the cuts remain non-intersecting. In particular from 
now on we will choose the cuts so as to lie on the triangulation. This implies a restriction on 
the class of triangulations considered as such a choice may not exist for small triangulations. 

Now we can state: 

Lemma 2. 

Zpn{^,J^')= [ lldh,m%) (4) 

where S"^ is a handlebody, dT,^ is its triangulated boundary which carries half integer labels on 
its edges and C labels the cuts. Choose a set of cuts Dj that lie on the triangulation. (((9S'^)J^.) is 
the spin network evaluation of the dual of the triangulation of the surface, with the links labelled 
by the half integer lengths of the edges they cross and a hi G SU(2) inserted on every link that 
crosses a cut 5Dj G dH^, i £ C. The spin network is evaluated in the planar projection of the 
boundary of the cut manifold. That is, with all crossings being due to the links crossing a cut. 

Proof: Cutting T,^ along the discs Dj reduces it to a 3-ball. The spin network evaluation is 
defined by taking the planar no crossing representation of the graph cut along the circles dOi 
and then connecting the identified open ends. If we choose a triangulation that triangulates 
each disc Dj without internal vertices and reexpress the resulting amplitude as a spin network 
evaluation, then the gluing of the faces corresponds to a projection onto invariant subspaces. 
Replace the projection onto the invariant subspace by a group integration and we get the lemma. 
□ 

Note that due to the intertwining property of the spin network the choice of Dj does not 
matter, it merely moves the hi insertion in the intertwiner around. 

2.2 Coherent triangles 

In order to have a clear geometric picture of the amplitude we will choose the intertwiners 
in the boundary state ^' by using coherent states ak{n,6) [16]. These are the highest weight 
eigenstates of the normalized Lie algebra elements, that is for = |cj* the Lie algebra generators 
and n G S*^, a coherent state a/c(n, ^) in the k representation satisfies: 

L.nak{n,0) = ikak{n,6) (5) 

The parameter 6 describes a choice of representative of the U(l) equivalence class of states 
that correspond to the same n. These states transform with a phase under the group elements 
generated by L.n and the label n transforms covariantly under the S0(3) action of SU(2). That 
is for g G SU(2) with corresponding S0(3) element g: 

gak{ri,9) = e"'^ak{gn,e) (6) 

For the asymptotic analysis three further properties will be crucial: 
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• The k representation can be constructed as the symmetric subspace of 2k copies of the 
fundamental representation. In this picture coherent states decompose into a tensor prod- 
uct of coherent states in the fundamental representation. Consequently the group action 
factorizes: 

2k 2k 

gak{n,9)=g(^a^{ii,e) = l^e'2ai{gn,e). (7) 

i=l 1=1 

• The modulus squared of the Hermitian inner product of coherent states is given by: 

\{ak{ni,di),ak{n2,e2))\^ = l-{l + ni.n2)] , (8) 

• Under the action of the standard antilinear structure on SU(2) (see [9]) the coherent state 
changes as: 

L.nJak{n,0) = —ikJak{n,0) (9) 

The antilinear map J is given by multiplication by the epsilon tensor in the spin k repre- 
sentation followed by complex conjugation. J commutes with SU(2) elements. 

Note that given a set of three edge labels ki there is a non zero intertwiner exactly if they 
satisfy the triangle inequalities. Therefore there is a set of rij, unique up to 0(3) such that 
Ylf=i ^i^i — 0- We can then choose our intertwiner in the boundary state ^' as 

L= [ dX(^{Xak{n^,e,)) (10) 

JSU(2) Y 

This state is clearly an SU(2) invariant state. As we noted that SU(2) acts covariantly as S0(3) 
on the labels iij this choice is only dependent on an unspecified phase as we left open which 
eigenstates of L.Ui we are using. In particular it does not depend on the remaining parity 
P = 0(3)/SO(3) as this acts on the plane of the triangle as an S0(3) element. 

Thus choosing normalized ak{n,0) compatible with the boundary spin labels fixes the inter- 
twiner states up to a parity choice and up to a phase. These two data will be fixed by considering 
the gluing of the boundary. 



2.3 Regge state 

Let y be a set of labels for the boundary faces. Then we label the boundary edges by pairs 
ab \ a,b €z V and call the set of such pairs E. Let (j)a '■ — be an orientation preserving 
map from the a-th triangle on dT,^ to the plane orthogonal to the north pole of 5^ (which we 
denote = (0,0,1)). We choose the orientation in AA"*- to be the one inherited from by 
taking M to be the outward surface normal. As the boundary of is orientable, we can define 
Hafe = Meab) where ab G E. 

The requirement that (pa be orientation preserving implies that the triangles with the edge 
vectors given by /cjiij all have the same orientation in AA"*-. In particular we can require them 
to have the same orientation as we have chosen for Af-^. In particular this implies that we can 
glue up any two triangles a,b with a common edge in AA-*- in an orientation preserving way. 

Thus there exists an element gab G S0(3) such that: 

M = Qab^f (11) 

Where nab is the edge vector of triangle a that gets glued to triangle b. As in [9] , gab is the Levi- 
Civita parallel translation from triangle a to triangle 6, according to the bases provided by 0a 
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and (pi,. Again, given a choice of spin structure for S'^, a choice of a spin frame for each triangle 
defines the SU(2) lift gab as the parallel translation of the spin connection in these frames. 

Next we will describe a canonical choice of phase for the boundary state ^. From ([6]), 
afc(— riafej^ab) is proportional to gabOik{'!^ba-,0ba)- Then we fix the relative phase of the coherent 
states on the boundary: 

Jak{nba,Oha) = Qabaki^abidab) (12) 

We call coherent states with the above relative state choice Regge states, and denote them |n, k) 
Their image under the antilinear structure is |— n. A;) = J|n, A;), and states in the fundamental 
representation are denoted |n). 

The total boundary state is then given by: 



*(fei,ni) = 




Due to the presence of the antilinear map in the definition of the relative phase the overall 
ambiguity not fixed by (I12p cancels in the overall state. At each triangle a we have a sign 
freedom as adding a sign contributes {^—\)'^'^b,ab£E^a.b _ ]^ ]-,y ^-j^g admissibility conditions on 
intertwiners. This shows that as in [9] the possible lifts of gab are defined by the spin structures 
on the boundary, and do not rely on the arbitrary spin frame covering chosen to define the lift. 

Finally note that inverting the orientation of A^-*- would have the same effect as turning the 
state into ^' = . 



2.4 The Amplitude 

We begin with i?'^. To evaluate the spin network defining our amplitude in terms of these 
coherent intertwiners we choose a particular diagrammatic representation of the planar graph. 
To obtain the spin network evaluation of this graph we then contract the intertwiners chosen 
using the epsilon inner product defined in terms of the Hermitian inner product by (a, (3) = 
{Ja\(3). Number the triangles in the graph from left to right. We then assume that the coherent 
intertwiners have been specified with respect to this planar representation of the graph as well. 
Then we have no crossings in the diagram and we can now explicitly write the contraction of 
coherent intertwiners as: 

Zpr{^, B^) = j W dXa n (^fe |nbc, he) , X, |n,b, he)) 

aGV be£E 



I 

[HdXaH (-ribel Xlx, M^'"^ (14) 



Y[ dXa Y\ (-nfoc, he\ XlXc Iricfe, hb) 
a£V beeE 



aGV beeE 

Where we have written \ncb) for \ncb, 

For general manifolds we need to make sure that, after we have chosen the circles, the dual 
edges crossing a circle all have the same orientation relative to the circle. This can be done 
by using a planar representation that has all the glued discs strictly left or right of each other. 
Call E the set of edges not crossing circles and Ej the set of edges crossing circle j £ C. The 
amplitude is then given by 



Zpni^, S3) = (-l)x JUdXall dh, II (-ribj Xlx, |n,b)2'=^= J] (-n^,] X%X, |n 



edi 



aeV jeC bc&E l€C deeEi 

(15) 
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Where (—1)'^ is a sign factor incurred in the spin network evaluation when connecting up the 
glued edges in the spin network evaluation. This can then be written as 

iev jec 

with the action given by 

S=Y, 2U In {nab\ JXlXb \nba) + Yj Yl ^^^^ ("'^^l J^d^iXe M- (17) 

abeE l&C deeEi 

Note that the ambiguity in the logarithm of a complex number does not affect the amplitude. 



2.5 Symmetries of the action 

The action ()17|) has the following symmetries (up to 27ri) 

• Continuous. A global rotation Y G SU(2) acting on each Xa and hi as Xa YXa and 
hi — > YhiY~^. This represents a rigid motion of the whole manifold. 

• Discrete. At each triangle a the transformation Xa — > CaXa with = ±1 leaves a factor 
Ca '''"'^^ " . As the admissibility conditions are satisfied on each triangle, this factor 
equals one. Similarly we have an arbitrary sign on hi as the edges on which hi act 
satisfy the admissibility condition for intertwiners. 

This latter symmetry will be used to compensate for the ambiguity of the lifts of S0(3) to SU(2). 
We can now state the theorem on the asymptotic formula. 



2.6 Relation to the standard intertwiner phase choice 

The standard choice of phase for an intertwiner, defined by chromatic evaluation |T7], gives real 
numbers for a spin network evaluation. We will now show that with the Regge phase choice the 
amplitude is real so can only differ from the chromatic evaluation by ±1 and a normalisation 
factor. Note that since the Regge choice has all the n^fe orthogonal to e^, the rotation g"*'^''^''^ 
rotates ricb to —n^b and leaves invariant. Under this rotation, the coherent state In^b) will 
transform as 

e-'-'^-nncb) = e''^J|n,fe) (18) 

for some phase 4>. 

Consider a single term in the amplitude (1151) . and rewrite it inserting the identity: 



{-nb,\ Xlx, M = {-^bc\ e*™-'^(e-*™-'^X,t)(Xee*™-'^)e-^-^-'^ M 

= (n^el J^e-'^J^xlx,g,bJe''^J\^bc) 
= {Mbc\J^J^Xlx,J\ricb) 

= {-r^bc\xlx,M (19) 

where we have defined the transformation X^, = Xce^'^^'''"' , which can be absorbed on the group 
integration in (fTSl) and the fact that J \ncb) = \—ncb). We have used the Regge phase choice (fT2]) 
from going from the first to the second line and the fact that the SU(2) transformations are all 
in fact in the same U(l) subgroup (and hence commute). From going to the second to the third 
line we have noted that we are acting with opposite rotations on the same state. Hence we get 
that Zpr{^, E^) = Zpr{^,T,'^) which is thus real. 
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h G S0(3) 



Figure 3: A cut immersion for a particular boundary triangulation of a torus. The cutting circles 
are shown in bold and there is an /i E S0(3) that identifies them. 



3 Asymptotic formula 

We wish to study the semiclassical limit of the amplitude Zpi^(^', S^). In order to do this, we 
homogeneously rescale the spin labels by a factor A. The corresponding boundary state is 
given by *a = ^(AA;^, rii). 

Given a set = {iiab, kab}ay^b of boundary data we denote as J the set of cut immersions of 
the polyhedral surface (9S^ with edge lengths kab in up to rigid motion. 

A cut immersion i G 21 is an immersion of the manifold obtained from dT,^ by the trivializing 
cuts dDi, i € C, i.e. it is an immersion i{dE^ - {Uigc^Oj}) ^ M^. Furthermore, we require 
the existence of S0(3) elements that identify the two sides of the cut, i.e. hi S S0(3) such 
that hi{i{dOf)) = i{dB~), where dB~ and dBf are the elements of the boundary — 5Dj) 

created by the removal of dDi from dT^^ . 

Any two cut immersions of are defined to be equivalent, and can be obtained from each 
other if the cuts are related by a homotopy on the surface. Therefore, different choices of cuts 
Bj lead to equivalent cut immersions. An example of a cut immersion is given in Figure El 

Such an immersion is called rigid if every continuous deformation of it requires changing the 
edge lengths, and flexible otherwise. We denote the subset of rigid immersions 3r C J. Through 
every immersion in J passes at least one manifold (with dimension d) of immersions that can be 
continuously deformed into each other. We call these flexifolds and denote them f, we denote 
the set of flexifolds ^. We then define 5max to be the set of flexifolds in ^ of maximal dimension 
c^max- With this definition the rigid immersions are a special case of a flexifold with dimension 
d = 0. We assume from now on that the flexifolds f do not intersect. 

In the limit A — > oo we have the following theorem: 

Theorem 1. (Asymptotic formula) 

1. If 3 is not empty we have that: 



S(\V\ + \C\-l)-dmax / 



S(\V\ + \C\-l)-dmax 

S(\V\ + \C\)-d,n ax 

+0\{\) ' I (20) 



^Note that since Di is transversal to generators of _H"i(S''), its removal changes the connectivity of and 
creates two new boundaries, D~ and Hf. 
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The coefficient Nf, the dihedral angle Q'^^ and the phase (t)\^ are independent of A. The 
and the dihedral angle 0^^^ are evaluated on an arbitrary immersion i m f. It can he 
shown that these are independent of the cuts. Thus for any particular edge we can evaluate 
the dihedral angle by moving the cut away from it. \V\ is the number of triangles (or 
equivalently vertices in the set V) and \C\ is the number of cutting circles. is the 

dimension of the flexifolds f G S,nax, (^nd Nj now also contains an integral over the union 
of flexifolds in f. 

2. If no immersions in exist the amplitude is exponentially suppressed: 



ZpRi^x,^^) = o{\-^) Vn (21) 

Note that in the simple case where the boundary data only admits rigid immersions, ie if 
c^max = 0, then the sum becomes a sum over the rigid immersions \ and we have that: 



3(|Vi + |C|-l) / 

Zpr{^^,t?) = (^) ' ^iv^cos A^uel, + <Ai 



A 7 

iGdr \ ab&E 



ab 




{22) 



since dmax = 0. Since the immersions are now rigid, the coefficient A'';, the dihedral angles 
and the phase (/>Jjj are evaluated on the cut immersion i. 



4 Proof of the asymptotic formula 

We now prove the above theorem. We begin by describing the methods used to give the asymp- 
totic form of the amplitude, this will require finding the so-called stationary and critical points 
of the action. We can then interpret these points geometrically and give the asymptotic formula. 
Much of this section is similar to [9] but one dimension lower so the analysis will be as brief as 
possible. 



4.1 Stationary phase 

To find the asymptotic form of the amplitude /, we will use the complex stationary phase formula 
|18j . The details of this are recalled below. 

Let D be a closed manifold of dimension n, and let S and a be smooth, complex valued 
functions on D such that the real part ReS" < 0. Consider the function 

/(A)= / dxa(x)e^^(^). (23) 
Jd 

The Hessian of S is the n x n matrix of second derivatives of S denoted H. For now let us assume 
that the stationary points are isolated, which, by the Morse lemma is a condition equivalent to 
the statement that the Hessian is non-degenerate at the critical points; det H ^ 0. Such functions 
are called Morse functions. 

In the extended stationary phase, the key role is played by critical points, that is, points xq, 
which are not only stationary: 6S{xo) = but for which Re5(xo) = as well. So to compute 
the dominant terms in the asymptotics for large spins, we need to find the stationary points 
of the action S and restrict to those with zero real part. If S has no critical points then for 
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large parameter A the function / decreases faster that any power of A .In other words, for ah 
N >1: 

/(A) = o(A-^), (24) 

If there are isolated critical points, then each critical point contributes to the asymptotics of 
/ by a term of order A~"/^. For large A the asymptotic expansion of the integral yields for each 
critical point 

(r) \ n/2 -J 

At a critical point, the matrix —H has a positive definite real part, and the square root of 
the determinant of this matrix is the unique square root which is continuous on matrices with 
positive definite real part, and positive on real ones. If S admits several isolated critical points 
with non-degenerate Hessian, we obtain a sum of contributions of the form ()25l) from each of 
them. 

If there are also degenerate critical points, more care is needed to compute the asymptotics 
(see eg [E]). Let C := {y £ D \ 6S{y) = 0, Re<S'(y) = 0} denote the set of critical points. For the 
action S : SU(2)I^I"''I'^I ^ C we will show that C is the set of immersions 3. Note that we cannot 
a priori assume C to be a disjoint union of manifolds as the flexifolds in 3 can intersect. However, 
here we have restricted ourselves to this case so that the following generalized stationary phase 
theorem applies. 

For a smooth function S whose critical set C is a disjoint union of closed manifold^, each 
critical manifold Cxq of dimension p, labelled by some xq on the critical manifold, contributes 
the following to the asymptotic formula |19j : 

where (y) is the restriction of the matrix to the directions normal to Cx^ with respect to some 
Riemannian metric on the domain , and dar is the canonical measure induced on the critical 
submanifold by the same Riemannian measure on the domain space. This extends to the case 
where C is a manifold-with-boundary. 

4.1.1 Critical points 

As described above, we must find the points of the action (|17p such that ReS* = as these are 
the only points that contribute in the limit A oo. First, we introduce some more notation. 
The action of the elements Xi, on the coherent states will produce a new coherent state 

K,) = (27) 

We will denote the corresponding rotated three vectors by 

n^fe = Xanab (28) 

where Xa is the S0(3) element corresponding to Xa- 

We will first consider critical points for edges that are not on one of the cutting circles. Using 
([H]), we can see that the real part of the action is given by 

Re5= Y,kabln^il-K,-nl). (29) 

abGE 

^In the literature this is called a Morse-Bott function. A Morse function is the special case where the critical 
manifolds are zero-dimensional (so the Hessian at critical points is non-degenerate in every direction, i.e., has no 
kernel) . 
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This does not depend on the coherent state phases as it is real. Using this formula, we can see 
that ReS = when n^^^ = — n^^j for all ab, or explicitly in terms of S0(3) rotations 

Xallab = -Xb^ba- (30) 

The critical points for an edge that crosses a cutting circle i differ by the inclusion of the hi 

Xa^-ab = —hiXhUha- (31) 



4.1.2 Stationary points 

The stationary points are found by varying the action with respect to each of the group variables 
Xa- The variation of an SU(2) group variable and its inverse is 

5X = TX = -X~^T (32) 

for an arbitrary su(2) Lie algebra element T = ^iT^aj. The stationary points are given by 
5S = and lead to the following equation 

kab V,fe = (33) 

b: bjta 

where 

V,, = i-^ab\X-' ^ Xb\nb.) ^ ^^^^ 

{-nab\Xa Xb\nba) 

These equations can then be evaluated at the critical points to give 

kab nab = (35) 

b: b^a 

which is the closure constraint for an immersed triangle. 

The stationary phase condition for the hi variables is the same but in this case we obtain 

E <b = (36) 

ab&d 

Which is the closure condition for edges on the circle i immersed in . Note that unlike the 
closure condition for the triangle, this relation involves the n^^^ as each edge belongs to a different 
triangle. 

If the critical points are not isolated but form a manifold of critical points, then we denote 
this manifold by 

Cx = {(Xi, ...,X\y\,hu h\c\)^ SU(2)I^I+I^I : 6S = 0,Re(5) = 0} (37) 



4.2 Geometry 

We will now describe how the critical/stationary points can be given a geometric interpretation. 
We first consider the case where the critical points are isolated. 

Theorem 2. (Geometry) Given a set of boundary data B satisfying the closure constraint on 
each triangle, the solutions Xa,hi to the critical and stationary point equations (j30p . (j3ip and 
(I36p correspond to the oriented geometric immersions \ of a geometric triangulated 2-manifold 
with boundary (Jiec ^^/^ ^ obtained by a suitable cutting of the boundary manifold. This 

immersion is subject to the constraint that a set of hi E S0(3) exist that map the immersion of 
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dnf to the parity flipped P((9D^), that is, the immersion of dB)f is congruent and oppositely 
oriented to the immersion of dO~ . 

The geometric vectors of the immersion are given by 

Va6(i) = kabXallab- 

Its orientation the one induced by the vectors on each face. 

Conversely, an immersed surface i determines a set of kab,nab and a set o/SO(3) elements 

Proof: Start somewhere on the surface of that is not in a circle. Since has connected 
boundary, the entire surface is now contractible to this point if cut along the circles. Take the 
triangle Aq you are on as embedded in M'^ and rotate it according to the stationary point 
equations to XaAa- Embed the next triangle and rotate it, according to the stationary point 
equations it's edges are now antiparallel to already immersed edges. As they are geometrically 
glued a translation exists that identifies all its edges with already immersed ones. Thus iteratively 
the whole immersed surface can be built up and the closure conditions on the cuts now imply 
that the surface closes up. Finally the stationarity equations on the circles imply that the Phi 
identifies the circles where we cut the surface. 

Conversely given an oriented immersion with the right edge lengths we can choose a set of 
edge vectors compatible with the orientation on the surface. On each triangle there are two 
linearly independent edge vectors. The map from these to the corresponding boundary elements 
defines a rotation in S0(3). On the boundary circles we explicitly have S0(3) elements. The 
complete set of these solves the critical point equations. □ 

If there is a manifold of dimension d > of critical points then Theorem [2] holds for each 
critical point in Cx- Since these critical points lie on a manifold, there is a continuous deformation 
of the immersed surface that does not change the edge lengths. Hence these critical points 
reconstruct flexible immersions and we arrive at the flexifolds f described in section [31 We will 
now label the critical manifolds by Cf, where f is the flexifold that it describes. 

Geometrical interpretation We will now describe in some more detail the geometric struc- 
ture of these surfaces. Given an oriented surface in M"^ the standard orientation automatically 
gives us a consistent set of normals Ma- By our choice of boundary data we have automatically 
ensured that these are given simply by 

Ma = XaM. 

We can define the dihedral rotation for an oriented surface unambiguously as the rotation Dab £ 
S0(3) around the geometric edge vector ^ab{^) that takes Ma to Mb, that is 

Mb = DabMa 

and 

Vab{\) = Dab^abii)- 

the lift of this rotation can thus be written as 

Dab = expfei^^.L) 

= exp(Gl,<,.L) (38) 

where we require — vr < GJ^^ < vr. We then call OJ^^, the dihedral angle. As Vabi^) = — Vfea(i) this 
definition clearly implies GJ^^ = G^^. If we have a surface defining a convex subspace of this 
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definition reduces to the usual definition. In particular the consistent choice of orientation then 
ensures that we have < < vr and cos (0Jj^) = Ma - Mb for outward facing normals. 

On the boundaries of the surface we can also define an analogue of the dihedral rotation by 
requiring 

hiMb = DabMa 

and 

Geometrically this makes sense as it corresponds to the angle obtained by gluing on the two 
identified boundaries. 



Commuting diagram To fully connect the geometry of these surfaces to the boundary we 
can relate the dihedral rotation to the gluing of the boundary gab defined above. 

Consider the following diagram which applies to two adjacent triangles that are not on a 
cutting circle: 

Xa 



9ab 




(39) 



Here ta is the boundary triangle at AA-*- with edge vectors given by kab^-ab and Ta is the 
triangle rotated according to its location in the surface, which according to the reconstruction 
theorem [2] has edge vectors given by Vafe(i) = kabXanab- The SO (3) action of this diagram 
immediately commutes, as can be seen by acting on rvab and N . As an SU(2) diagram this 
equation defines a sign (—1)'''''' that makes it commute. The discrete sign symmetry Xa ^aXa 
of the action can be seen as acting on this sign by (—1)'''''' 1)'^"''. 

Now, for two triangles whose common edge is on a cutting circle i, in the same way we have 
a commuting diagram as 



9ab 




(40) 



(_l)-ai,D,i 



and additionally have {—\Y^*> — > eaebej(— l)'^"*. 

We now show that the dihedral rotation is unchanged by moving the cut. Suppose that 
af) G i G C, in the sense that the edge yrab{^) crosses the cutting circle labelled i. Now choose 
a different, yet homotopic, cutting circle (?Bj ~ (9Dj, such that now ba ^ i, but cb G i. This 
corresponds to sliding the cut from one edge of Tb to another, and as any other change in the 
cut, it will modify wether labels a and b satisfy (I30p or ()3ip . By our orientation convention, we 
now have 

Xa'O-ab = —Xb^-ba , — hiXcTClcb = Xblilbc (41) 

By just looking at the first set of equations above, and comparing it to Xanab = —hiXbiiba we 
have that hiXbUba = Xbiiba- This shows that hiXb and Xb are equal up to a phase, and by the 
reconstruction theorem in fact hiXb = Xb- Hence, comparing the diagrams (I39p and ()40p we get 
that moving the cut does not affect the dihedral angle as here defined. 
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Parity Given an oriented surface immersed in M"^ all surfaces related to it by 0(3) clearly 
have the same boundary data. Those related by SO (3) also have the same dihedral angle as 
defined above. However if we act by parity P : n ^ — n we switch the dihedral angle. This is 
because the two equations defining Dab are invariant under parity, and the dihedral rotation is 
unchanged. Thus by the definition of the dihedral angle we have 

Dab = eKp (Q[^(-^{Pa)).L) 

V IVabI / 

and so ©S,^ = 

4.3 The Regge action 

Now note that given a geometric surface and associated solution Xa (i) , we can obtain the solution 
Xa{Pi) corresponding to the parity flipped surface explicitly. That is, by equation ([38|) we have 
that 

{Xa{x))-'DabXa{x) = exp n.^.L) (42) 

Acting with X~^ by the left of the commuting diagram equations, using the notation Xab{^) = 
{Xa{\))~^Xb{V) if not on a circle and ^^^(i) = {Xa{V))~^ hiXb{V) if on, we get 

Xabi^Qab = i-lY'^' {Xa{x))-^DabXa{i) = (-l)'^"" exp {@ab {^ab-L). (43) 

where we have used (j42p . 

It is now straightforward to evaluate the matrix elements in the amplitude. These are of the 
form {-Wab\Xab\Ta.ba)- Using the gluing condition this becomes {-Wab\Xabgab\-Tn.ab) ■ Finally, 
by ()43p this is just (— l)'^''''e2®"*'. Thus we have overall that: 

{-TYab\ Xab Ha) = (-l)'^-^et®ii> (44) 

Spin Structure Now we will fix the ambiguity in signs emerging from the spin lifts of the 
dihedral angle by exploring the discrete sign symmetry in hi and Xa- Recall that the discrete 
sign freedom of the action Xa — > ^a^a emerged from a different choice of spin frame for each 
triangle. Now, we show that the discrete sign symmetry related to the cuts hi — > ej/ij corresponds 
to different choices of spin structures for the manifold . Then, using the fact that {—lY'^^Dab 
is a gauge transform of the connection gab we can fix the symmetries by adjusting the spin frames 
and the spin structure such that (—1)''"'' = 1. Thus we will show that: 

Lemma 3. The signs (—1)'^"'' arising from the spin lift on each face not on the cut obey (—1)'^"'' = 
i^ab = KaK,b for some Ka = ±1. The signs for a face on the cut, i.e. ab £ i £ C obey {—1)'^"'' = 
Kab = KaKbKi whcre Ki paramctrizcs the spin structures ofTj^. 

Proof. First of all, by ()4.2p a lift of the dihedral rotations, HabDab-, are just a gauge trans- 
formation of the gab- Now recall that gab G S0(3) are parallel translations on the boundary 
triangles according to the Levi-Civita connection of the associated metric, with gab being the 
parallel translation of the respective spin connection (a lift of gab to SU(2)). 

But when the geometry around a vertex is continuously deformed to the flat geometry, the 
gab holonomy of a trivial cycle around said vertex has to go to the identity rotation, as opposed 
to a 27r rotation. This implies that for the holonomy around a vertex through triangles a, b and 
c (which of course consists of a trivial cycle), we have 

l^cal^bc^ab^-^ caDbcD ab — f^caf^bc^ab — 1 
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which imphes that locally we must have = i^ai^b- The problem now is that if there are 
non-trivial cycles, i.e. g' 7^ 0, we may not be able to extend this globally, i.e. Kab may not be 
globally pure gauge. 

In other words, for trivial cycles the lift of the holonomy given by the KabDab is fixed to 
be the same as that given by Dab- But not so for the holonomy of a non-trivial cycle; there 
exist inequivalent spin structures on a manifold. These have a one-to-one correspondence with 
the elements of -ffi(S^,Z2), and so are 2^ in number. Hence for a non-trivial cycle, dual to the 
sequence of triangles A^o • • • Aq^Aqq crossing the circle i € C, we have 

'^a„ao'^^a„_ia„ " " " ^aoai-Da„ao-^a„_ia„ " " " Da^ai = l'^iDanaoDa,^_^an ' ' ' ^aoai (45) 

where a Kj is introduced whenever there is an implicit choice of spin structure; i.e. it parametrizes 
the different spin structures associated with the cut. 

We reconcile this case with the g = one by keeping the form Kab = i^a>^b for all the edges 
ab that do not lie on a circle, i.e. ah ^ i for any i £ C. Then by (I45p immediately we must have 
for ab G i, Kab = i^ii^af^b- Since our chosen basis for HiiTi^ ^^,2) generates all cycles, we can see 
that this form of Kab has all the right properties demanded by our equations and accounts for 
the different spin structures. 

Therefore, taking advantage of the discrete sign symmetry, we can choose the spin structure to 
be compatible with the one chosen for the lift of gab and thus we will have (—1)'^'''' — > eieaeb{—'iY'^^ 
makes {-ly^^ = 1. □. 

Finally, we obtain that the action evaluated at the critical points is the Regge action for the 
immersed surface i. 

S=Y. kabOU (46) 

abeE 

For the flexible immersions, the action is the same for all point^ on the critical manifold so we 
evaluate it on an arbitrary immersion in the flexifold. 
This concludes the derivation of the Regge action. 

4.4 Hessian 

The stationary phase formula requires us to calculate the Hessian of the action S to determine 
the weights with which the stationary points contribute to the action. This will be a 3{\V\ + 
\C\) X 3(|F| + |C|) matrix defined by 

H=( ^-^^ ""J'). (47) 



Where 



'"-'-n^S) '''' 



^This is actually a particular case of the "strong bellows conjecture" that was shown in 
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The global SU(2) symmetry of the action implies that there is a redundant integration in /. 
This will cause the determinant of the Hessian to be zero unless it is gauge fixed. To solve this, 
we make the change of variables Xa X^Xa for some b S {l,...,|y|}. This has the effect of 
removing the Xfj variables and the integral gives a volume of SU(2) which can be normalised to 
one as it is compact. The remaining Hessian is now a 3(|1^| + \C\ — 1) x 3{\V\ + \C\ — 1) matrix. 
The submatrix Hxx at the critical points, is given by 



dXidXi 



1 



55=0,ReS'=0 b^c,bcGE 

for the diagonal terms. The off diagonal part is 



cb^ cb 



(52) 



dxi^dxl 



1 



55=0 
RcS=0 



c s{e)^d t{e)+^d s{e)^c t(e) 



ijk k 



s(e)t(e) 



^' s{e)t{e)n''i(^e)t{e) 



(53) 

So one can see that only the off-diagonal elements that represent two neighbouring triangles are 
non zero. The (H^h) submatrix will be diagonal since each term in the action only contains one 
hp term (ie, each dual edge only crosses one cut.) 



dh^dhi^ 



1 



^ J2 hbis'^-n'iyi, 



(54) 



5S=0,RcS=0 b^c,bc&Cp 

The mixed terms Hxh^Hxh will be non zero only for triangles with an edge on the cut 



dXidhi 



1 



55=0 
Re5=0 



ab&Eq 
c=a,b 



^£ ^ab 



n 13 ' 

ab'f^ ab , 



(55) 



Note that 

{Xa{x))-^DabXa{x)=ei^V{-Kb M<))-L)~^ = {{Xa{Pi))-^ DabXa{Pi))~' 



(56) 



where we used that Dab = exp {Q^^, (— |^^^|^^).L) on the second equality. By (fl3l) we then have 

{Xab{i)gab) = iXabiP^)gab)~^, SO if we replace the ^^^(i) in (-n^fel ^^^(i) \nba) with the parity 
related one we now obtain the complex conjugate: 

{-nab\Xab{i)gab\-T^ab) = {-T^ab\ {Xab{Pi)gab)~'^ \-T^ab) 



{-Uabl XabiPi)gab\-nab) 
{-nab\ Xab{Pi) \nba) 



(57) 



Thus we can see that the action of parity on the Hessian matrix will also result in complex 
conjugation when evaluated at the critical points. 



4.5 Proof of the formula 

We can now evaluate the stationary phase approximation to the amplitude Zpji{^\, S'^) defined 
in p5|) . We begin by fixing the symmetries of the action. This can be achieved by taking an 
arbitrary vertex and dropping the group integration associated to it. As shown in section [4.1.11 
the critical point equations are the equations for the immersion of a polyhedral surface with 
the geometry specified in the boundary data. If a particular immersion is rigid no infinitesimal 
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deformation taking it to another such immersion exists and therefore it is an isolated critical 
point of the amplitude. 

For the isolated critical points in 3r we can explicitly evaluate the stationary phase approx- 
imation. Having fixed one group integration we are left with a 3(|y| + \C\ — 1) dimensional 
integration. The overall scaling of these points is thus ^^^'^'^l^l"*"''-'! ^)/'^ ^ Further we obtain a 
set of 2^^l^l"'"l'-'l~^) critical points for each immersion from the spin lift of each SU(2). Finally 
the derivatives in the Hessian as defined above are taken with respect to a parametrization of 
SU(2) with volume (47r)^, so we need to rescale by this factor. Using equation (I44p and lemma 
[3] the amplitude itself evaluates to the Regge action of the cut immersion: 

In (-n^fel Xab Irifoa)^^"' = iKb^U- 

Taking all these factors together and using the fact that we know parity related immersions to 
be the complex conjugate of each other we can approximate the contributions of the isolated 
critical points to the partition function as: 

3(|V| + |C|-1) 3(|y| + |C|-l) / s 

where is the dihedral angle of the edge ah in the cut immersion i G J^- Since the Hessian 
matrix changes to its complex conjugate with parity, we can absorb the phase of the determinant 
into the exponentials and combine the terms 




(59) 

If there are any flexible immersions of the boundary data then there will be a manifold of 
critical points. Since the critical points extremize the action, it must have the same value on 
every point of the critical manifold. The Hessian therefore has zero modes along the directions of 
the flexifold and we must treat the integral as having further symmetries in the neighbourhood of 
the flexifold. Factoring out these changes the scaling of the contribution of these critical points 
by A''™'"'/^, where d^^^ is the dimension of the flexifold. Therefore these immersions dominate 
the rigid immersions if they exist, their contribution is given by: 

3(|y| + |C|-l)-dmax , _ , 3(|y| + |C|-l)-dn .,^ 



Lf exp ( iX kab@ii, I + exp ( -iX ^ kab&lh 

;5max V ab£E / \ ab&E / 




3(|V| + |C|-l)-dmax | ^ ' 



(60) 
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©^^ is the dihedral angle of the edge ab of a particular cut immersion i in the flexifold f. As the 
action is constant along the flexifold it does not matter where we evaluate it. is given by 

Lf= [ dac,{y)^^= (61) 



f 

where HJ is the Hessian matrix for the transverse directions which we can not give a general 
formula for. Combining the exponentials into cosines as above we obtain part one of the main 
theorem. 

Finally, if no immersions of the boundary data exist then there are no solutions to the critical 
point equations and the stationary phase formula gives that the amplitude is suppressed. □ 



5 Example: The Tetrahedron 

Here we apply the above results to the well known case of the asymptotics of the amplitude 
for a single tetrahedron which, with an appropriate choice of normalisation for the boundary 
intertwiners, will correspond to the 6j symbol. This is a special case of theorem [T] so the proof is 
the same as above. In particular, the critical and stationary point equations are the same and 
the action evaluated at these points reduces to the Regge action for a tetrahedron. Since the 
asymptotic formula for the tetrahedron is already known, we must verify that our formula agrees 
with this result. This also provides further evidence that the asymptotic formula for the 4d case 
derived using the same methods in [9] is correct. We begin by noting that, up to parity, the 
boundary data of a tetrahedron has only one immersion so the sum in the asymptotic formula 
disappears. 

We will also require an explicit formula for the Hessian matrix in order to perform the 
numerical calculations. For the tetrahedron, the Hessian is the 12 x 12 matrix defined by 



dXidXi 



HZ=\^:7-^\- (62) 



The global SU(2) symmetry of the action implies that there is a redundant integration in /. 
This will cause the determinant of the Hessian to be zero unless it is gauge fixed. To solve this, 
we make the change of variables Xa — > Xi^Xa for a = 1,2, 3. This has the effect of removing the 
X4 variables and the integral gives a volume of SU(2) which can be normalised to one as it is 
compact. The remaining Hessian is now a 9 x 9 matrix which, at the critical points, is given by 



2 

<55=0,ReS=0 b^c 



\dXidXi 

for the diagonal terms. The off diagonal part is 



\ E - n'l,n'i,) (63) 



dXidX^ 



h,, (6^^ - - n'l,n'i,) . (64) 



(55=0,Rc5=0 

5.1 Normalisation and scaling behaviour 

We can now compare our theorem with the Ponzano- Regge asymptotic formula for the 6j symbol. 
The Ponzano-Regge formula is 



\ki2 Xki3 A/ci4 
>^k23 A/c24 AA:34 j ^127rVol 



L=cos (^E(AU + ^)e., + ^^ 



(65) 
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where Vol is the volume of a geometric tetrahedron with edge lengths Xkab + 2 ^^^^ @ab are 
the dihedral angles. Note that the formula scales as A~^/^ due to the volume term. Currently, 
our formula for the tetrahedron contains the Regge action but the amplitude, phase term and 
scaling do not obviously agree with (j65|) . We will first consider the intertwiner normalisation, 
which will be necessary to obtain the correct scaling behaviour and some numerical factors, and 
then evaluate the Hessian numerically to check the agreement of the remaining terms. The main 
drawback of the coherent state approach occurs here as it is very difficult to obtain an analytic 
formula for the determinant of the Hessian matrix. 



5.1.1 Intertwiner normalisation 

For the 6j symbol, the three valent intertwiners are normalised by dividing by the square root of 
the theta spin network. The coherent intertwiners that we replaced these with, however, are not 
normalised. The normalisation of these intertwiners for the coherent tetrahedron was studied in 
|13] so we will briefly summarise the results for the case of the coherent triangle. 

The normalisation of the coherent intertwiner is given in terms of the three edge vectors of 
the triangle ni, n2, by the Hermitian inner product 



fA{ni,ki) = / dxY\{T^i,h\X\ni,ki 

JSV(2) TLi 



SU{2) 

dX exp 5a (66) 



/SU(2) 

where 

3 

SA = Y,2kM^^\X\^i)■ (67) 

This integral can be calculated exactly using [I6l [TT], the result being 

^ (1 - ni.n2)P(l - ni.n2)''(l - ni.n2)''(p + q)\{q + r)l{p + r)! 
•'^ ~ 2P+i+^{p + q + r + iy.p\q\r\ ' ^ ^ 

Where p = ki + k2 — k^, q = k2 + k^ — ki and r = ki + k^ — k2- 

The asymptotics of this intertwiner normalisation can also be found using stationary phase 
|13j . The stationarity of the action 5a gives the closure condition and the action evaluated on 
the critical points ±7 gives zero. 



3/2 

/A(ni, Xki) 



27ry/^ 2 1 



A J (47r)2 Vdet Ha 

^ (69) 



The additional factor 2 comes from the fact that both / and —I are critical points that give the 
same contribution to the action. Ha is the Hessian matrix of the action which is given by 



H 



^ dX^dXi 

= \Y.h{5^^ - n^.n^) (70) 
I 

We can now normalise our formula such that it agrees with the standard normalisation by 
dividing by a factor {fAaY^"^ fo'^ each triangle a. 
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Figure 4: Comparison of the 6j symbol (dots), the PR formula (red line) and equation (j7ip 
(blue dashed line) against the scaling A. The scaling factor A~^/^ has been removed to make the 
comparison easier at low spins. 



5.1.2 Numerical calculations 

With the intertwiner normalisations included in the asymptotic formula, we obtain 
Xki2 Xki3 A/ci4 ) Zpn{'^x,(^) 



Xk23 \k2i A/C34 j n;=i\/7A 



X cos 



[ XkabQab - ^Arg(det H) j (71) 

\a<b / 



Note that we have the correct scaling behaviour once the additional scaling factors from the 
intertwiners are included. The normalisation terms are real so do not contribute any additional 
phase. 

The formula for the equilateral tetrahedron with both the exact and approximate intertwiner 
normalisation was compared to the 6j symbol and the Ponzano-Regge asymptotic formula using 
Mathematica in Figure [H We see that our formula differs from the Ponzano Regge formula 
for low spins. The only point where our formula differs from Ponzano Regge is in the fact 
that the Ponzano Regge asymptotics are given in terms of the dihedral angles and volume of 
the tetrahedron with edge lengths Xkab + ^. Therefore the dihedral angles and volume change 
nontrivially with A. A stationary phase approximation extracts only the scaling behaviour with 
respect to lambda in the asymptotic regime and cannot register this type of low spin behaviour. 
This agrees as well as the PR formula for larger spins, however the agreement for very low values 
is not as good - Figure [H 
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Figure 5: A net showing a set of boundary data that reconstructs Steffen's flexible polyhedron. 

6 Example: Steffen's flexible polyhedron 

Here we discuss an example for which the second part of theorem[T]is relevant, that is we describe 
a set of boundary data that admits a flexible immersion. This particular example is taken from 
a flexible polyhedron with half integer edge lengths consisting of fourteen boundary triangles 
which was found by K. Steffen [20]. A net for constructing this polyhedron is given in Figure 
El and the corresponding spin network in Figure [6l Since Steffen's polyhedron admits a flex in 
one direction, we know that the flexifold is at least one dimensional. As a polyhedron, it is not 
allowed to self intersect but there may be other immersions with flexibility in more than one 
dimension. Applying the asymptotic formula with the same intertwiner normalisation as the 
tetrahedron in section O we would expect the scaling to be \~^'^/'^. 

7 Discussion and Conclusions 
7.1 Rigidity of cut immersions 

As discussed, the asymptotic formula produces a sum over all possible immersions of the bound- 
ary data in M^, including flexible ones. These flexible immersions scale differently and thus 
dominate the rigid immersions. The question of whether a particular polyhedron is rigid is a 
difficult long standing problem in mathematics. A classic result is that convex polyhedra are in 
fact rigid, however this does not extend to immersions and non convex polyhedra where counter 
examples, like Steffen's polyhedron discussed above, are known. 

If the boundary data is topologically S"^ then a theorem by Steinitz [21] applies that states 
that any simplicial complex with underlying space homeomorphic to a 2-sphere admits a sim- 
plexwise linear embedding into whose image is strictly convex. This embedding will indeed 
be rigid and we can conclude that for the ball 3^ will always be non-empty. 

To our knowledge the only more general results on rigidity of immersions are those giving 
conditions on bar frameworks, that is a graph immersed in W^, to be generically rigid. A bar 
framework is considered generic if the coordinates of the vertices are algebraically independent 
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Figure 6: The spin network corresponding to Steffen's flexible polyhedron. 



over the rationals, that is, there is no polynomial with rational coefficients that has these coor- 
dinates as roots. A graph is generically rigid if all its generic frameworks are. A set of sufficient 
and neccessary conditions for a graph to be generically rigid are known \^22\ [23] • Unfortunately 
this does of course not cover our case with half-integer edge lengths. 

Concerning the rigidity of cut immersions, which can be seen as bar frameworks with addi- 
tional constraints, nothing is known. 



7.2 Surface immersions vs interior immersions 

With the asymptotic analysis performed above, we explicitly obtain a sum over immersions of the 
boundary data weighted by the cosine of the Regge action for the immersed surface. Previously, 
asymptotics of the Ponzano-Regge model for larger triangulations could only be considered by 
taking the product of the asymptotic formula for each 6j symbol. We now illustrate schematically 
that, in a simple example, that this is in fact equivalent to the asymptotic formula above. 

We will consider the case of two tetrahedra cJi , a2 glued along a common triangle A and use 
the boundary normalisation that agrees with the 6j symbol. The partition function then reads 

, , / Xki Xk2 Xh \ ( Xki Xk2 Xk^ \ 

ZpR{^x,<T,UAa2) = ^ ^j^^ A^^5 Xk,i\Xkr Xk^ Xk, j ^^^^ 

We write the asymptotic formula for the 6j in terms of the Regge action for a tetrahedron a 

a': Z ill} = '^(''Mi^s,) + .M-i^s,)) (73) 

Where several of the factors have been absorbed into the amplitude N. Asymptotically, this 
gives 

Zpr{^x,(Ti Ua CT2) = N1N2 (exp(iA(5^, + S^^)) + exp{-iX{S^, + 

+N1N2 {exp{iX{Sa^ - Sa^)) + exp(-iA(5'^j - S^^)) 
= A^iA''2exp(^ fcevr) (cos(A(5^,Ut<72) + cos(A(5^,UaP'T2)) (74) 

eCA 
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Figure 7: Two different possible immersions of tlie boundary data for two tetrahedra ai, o"2 glued 
on a common triangle A. 



Where Pa is the parity related tetrahedron and we have used the fact that the Regge action for 
two tetrahedra becomes 



Thus the formula gives a sum over the two different ways of immersing the boundary triangles in 
M^, see Figure [7]. If we now consider larger triangulations, possibly with interior vertices, edges 
or faces, then the immersed boundary surfaces we obtain are related to immersions of the entire 
triangulation by adding the immersions of the interior triangulation. Note that the internal 
vertices can be placed anywhere in which leads to the divergent factors that appear in the 
Ponzano-Regge state sum for these kind of triangulations. The two parts of the cosine appear 
because these immersions can be orientation reversing on particular tetrahedra. These will 
then contribute with negative dihedral angles to the interior action. Thus the integration over 
possible interior geometries can be split into parts corresponding to particular orientations on the 
interior tetrahedra. The terms of this sum will then correspond to the sum of terms obtained by 
multiplying out the sums in the interior cosines. On shell the interior action on these immersions 
is zero so the overall amplitude only registers the boundary contribution to the action. In 
particular interior configurations corresponding to different geometries and orientations on the 
interior, and thus to different terms in the cosine, contribute with the same phase. 

7.3 Boundary states 

We also note that it is possible to select a particular immersion in the sum by choosing a 
boundary state peaked around a particular set of dihedral angles, see for example [Ml (HI [25] , 
This boundary state also selects one overall orientation of the immersion which removes the 
parity related term in the asymptotic formula. For non-rigid immersions, the boundary state 
would also have the ability to select a particular configuration of the immersed surface which 
would stop these immersions dominating the integral. 

A possible problem with the boundary state is that while it selects an orientation for the 
boundary, it was not clear if the orientations of the interior tetrahedra behaved consistently. 
This was considered in [26] and our result also suggests that these do not cause a problem as 
the asymptotic formula does not register these orientations. 

7.4 Conclusions 

In this paper we addressed the problem of asymptotics of larger triangulations for the Ponzano- 
Regge model. By reformulating the partition function as a spin network on the boundary and 
then rewrote this amplitude using SU(2) coherent states. While this particular feature will not 
be available for non topological theories one could expect that in general boundary data will be 
not suppressed if it can be continued to a solution of the equations of motion on the interior. 




(75) 



ecA 
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The asymptotic formula contains a sum over immersions of the boundary data weighted by the 
cosine of the Regge action. Interestingly, Ponzano and Regge point out in [l] that the different 
possible immersions corresponding to 3-nj symbols should contribute to the asymptotics but did 
not obtain a concrete formula. The presented work opens up the possibility to do an exhaustive 
analysis of the classical limit of the Ponzano Regge model including correlation functions on the 
boundary. As such it can serve as a toy model and proof of concept for conceptual issues likely 
to arise in all background independent theories. 

Of further interest would be to consider in more detail how the asymptotics obtained here 
can be obtained from the "product of cosines" picture. In particular to shed light on the issue 
of causality and orientation in spin foam models. 

Interestingly, and unexpectedly, we found that spin networks contain some information about 
the rigidity properties of surfaces. The scaling properties of a spin network correspond to the 
maximum dimension of flexibility if the geometry to which it corresponds has any non-rigid 
immersions. 
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A Example of the Ponzano-Regge amplitude as a spin network 
on the boundary of the solid torus 

Here we give a simple example of Lemma E] on the solid torus T. A non-tardis (degenerate) 
triangulation of the solid torus with three tetrahedra is given by 




The two triangles with edges ki,k2, are identified. The Ponzano-Regge amplitude is given by 

We choose the cutting disc D to be the triangle ki,k2, k-^ and perform the cut that reduces T to 
the 3-ball. A net for constructing the triangulation on the boundary is given by 





kji 








kj^^^^^^ 




ki 


k2 


k^ 




k2 


ks 




ks 
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The Ponzano-Regge amplitude can be expressed as the following spin network evaluation on the 
boundary, with a group integral inserted on each of the dual edges that cross D. 

k2 



^pM^,T) = /gu(2) dh 




k3 



Expressing this spin network in terms of 6j symbols gives equation (|76|) . 

References 

[1] G. Ponzano and T. Regge, "Semiclassical limit of racah coefficients," in Spectroscopy and 
group theoretical methods in physics (F. Block, ed.), pp. 1-58, North Holland, 1968. 

[2] J. Roberts, "Classical 6j-symbols and the tetrahedron," GEOM.TOPOL., vol. 3, p. 21, 
1999. 

[3] R. Gurau, "The ponzano-regge asymptotic of the 6j symbol: an elementary proof," 2008. 

[4] S. Garoufalidis and R. van der Veen, "Asymptotics of classical spin networks," 2009, 
0902.3113. 

[5] J. W. Barrett and C. M. Steele, "Asymptotics of relativistic spin networks," Classical and 
Quantum Gravity, vol. 20, p. 1341, 2003. 

[6] L. Freidel and D. Louapre, "Asymptotics of 6j and lOj symbols," Class. Quant. Crav., 
vol. 20, pp. 1267-1294, 2003, hep-th/0209134. 

[7] M. Dupuis and E. R. Livine, "Pushing Further the Asymptotics of the 6j-symbol," Phys. 
Rev., vol. D80, p. 024035, 2009, 0905.4188. 

[8] J. W. Barrett and I. Naish-Guzman, "The Ponzano-Regge model," Class. Quant. Grav., 
vol. 26, p. 155014, 2009, 0803.3319. 

[9] J. W. Barrett, R. J. Dowdall, W. J. Fairbairn, H. Gomes, and F. Hellmann, "Asymptotic 
analysis of the EPRL four-simplex amplitude," 2009, 0902.1170. 

[10] J. W. Barrett, R. J. Dowdall, W. J. Fairbairn, F. Hellmann, and R. Pereira, "Lorentzian 
spin foam amplitudes: graphical calculus and asymptotics," 2009, 0907.2440. 

[11] L. Freidel and K. Krasnov, "A New Spin Foam Model for 4d Gravity," Class. Quant. Grav., 
vol. 25, p. 125018, 2008, 0708.1595. 

[12] J. Engle, E. Livine, R. Pereira, and C. Rovelli, "LQG vertex with finite Immirzi parameter," 
Nucl. Phys., vol. B799, pp. 136-149, 2008, 0711.0146. 



26 



[13] E. R. Livine and S. Speziale, "A new spinfoam vertex for quantum gravity," Physical Review 
D, vol. 76, p. 084028, 2007. 

[14] S. Speziale, "Towards the graviton from spinfoams: The 3d toy model," JHEP, vol. 05, 
p. 039, 2006, gr-qc/0512102. 

[15] J. Moussouris, Quantum models of space-time based on recoupling theory. PhD thesis, St. 
Cross College, Oxford, 1983. 

[16] A. Perelomov, Generalized coherent states and their applications. Springer- Verlag, 1986. 

[17] L. H. Kauffman and S. L. Lins, Temperley-Lieb recoupling theory and invariants of 3- 
manifolds. Princeton University Press, 1994. 

[18] L. Hormander, The analysis of linear partial differential operators I. Springer- Verlag, 1983. 

[19] P. Ramacher, "Singular equivariant asymptotics and the moment map i," 2009, 0902.1248vl. 

[20] K. Steffen, "A symmetric flexible Connelly sphere with only nine vertices." See: 
http://demonstrations.wolfram.com/SteffensFlexiblePolyhedron/ for Mathematica code. 

[21] H. Cluck, "Almost all simply connected surfaces are rigid.," in Lecture notes in mathematics 
438: Geometric Topology, Springer- Verlag, 1975. 

[22] R. Connelly, "Generic global rigidity," Discrete Gomput. Geom., vol. 33, no. 4, pp. 549-563, 
2005. 

[23] S. J. Cortler, A. D. Healy, and D. P. Thurston, "Characterizing generic global rigidity," 
2007, 0710.0926. 

[24] C. Rovelli, "Graviton propagator from background-independent quantum gravity," Phys. 
Rev. Lett, vol. 97, p. 151301, 2006, gr-qc/0508124. 

[25] E. Bianchi, E. Magliaro, and C. Perini, "LQG propagator from the new spin foams," 2009, 
0905.4082. 

[26] E. Bianchi and A. Satz, "Semiclassical regime of Regge calculus and spin foams," Nucl. 
Phys., vol. B808, pp. 546-568, 2009, 0808.1107. 



27 



